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ABSTRACT 

Let  P  and  Q  be  two  disjoint  simple  polygons  having  m  and 
n  sides  respectively.  We  present  an  algorithm  which  detennines 
whether  Q  can  be  moved  by  a  sequence  of  translations  to  a 
position  sufficiently  far  from  P  without  colliding  with  P,  and 
which  produces  such  a  motion  if  it  exists.  Our  algorithm  runs  in 
time  0(mn  a(mn)  log  m  log  n)  where  a(^)  is  the  extremely 
slowly  growing  inverse  Ackermann's  function.  Since  in  the 
worst  case  n(m/i)  translations  may  be  necessary  to  separate  Q 
from  P,  our  algorithm  is  close  to  optimal. 

1.  Introduction 

In  this  paper  we  develop  an  algorithm  for  the  problem  stated  in  the 
abstract.  This  problem  generalizes  previous  research  on  translational 
separability  of  planar  objects  (see  Toussaint,  1985,  for  a  survey  of  this 
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research).  In  most  of  the  previous  work  on  planar  separability  the  goal  was 
to  separate  the  given  objects  by  translating  them  one  at  a  time  in  some  single 
fixed  direction.  In  the  case  of  two  simple  polygons  P  and  Q  as  above, 
Toussaint  (1986)  has  given  an  0{m  +  n)  algorithm  for  determining  whether  P 
and  Q  can  be  separated  by  a  single  translation  of  one  of  them  (see  also  Sack 
and  Toussaint,  1985).  :'o;nv/  not^"'.).  '. 

In  this  paper  we  study  the  prcblenrnof;  sei5ar^tin[j  Q  from  P  under  an 
arbitrary  sequence  of  translations  (ijcv-  Fig.  1.1  for  an  illustration  of  this 
problem).  This  version  of  the  problem  is  of  intermediate  complexity  between 
the  simpler  single-translation  separability  problems  mentioned  above  and  the 
more  difficult  problem  of  separating  Q  from  P  by  any  collision-free  motion, 
involving  both  translations  and  rotations.  In  fact,  the  problem  that  we  study 
is  a  special  instance  of  the  motion  planning  problem  which  seeks  a  purely 
translational  collision-free  motion  of  a  polygonal  object  (Q)  amidst  a 
collection  of  polygonal  obstacles.  Our  case  is  special  because  the  obstacles 
consist  of  a  single  polygon  P.  Also,  the  destination  position  of  Q  is  fixed 
(some  placement  sufficiently  far  from  P);  however,  our  techniques  can  also 
handle  the  case  of  an  arbitrary  destination  of  Q . 

Several  related  motion  planning  problems  have  been  recently  studied.  If 
the  moving  object  Q  is  a  convex  polygon  (having  m  sides)  and  P  is  an 
arbitrary  collection  of  polygonal  obstacles  (having  n  sides  altogether)  then 
one  can  plan  a  purely  translational  motion  of  Q  between  any  two  placements 
(when  such  a  motion  exists)  in  time  0{mn  log  mn)  (cf .  Leven  and  Sharir, 
1985b,  Kedem  and  Sharir,  1985,  Kedem  et  al.,  1986,  Chew  and  Drysdale, 
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1985,  Fortune,  1985).  In  fact  within  this  time  bound  one  can  calculate  a 
discrete  representation  of  the  entire  space  of  free  placements  of  Q  (all  having 
the  same  given  orientation). 

Another  related  problem  is  the  polygon  containment  problem,  in  which, 
given  two  polygons  P  and  Q,  we  wish  to  determine  whether  Q  can  be 
transformed  into  a  polygon  which  is  entirely  contained  inside  P.  Chazelle 
(1983)  has  shown;  tbalijifii^is  dm^'ea  and  only  translations  of  Q  are  allowed, 
then  the  existence  of  such  a  placement  can  be  determined  in  0{m  +  n)  time. 
Related  work  on  polygon  containment  by  translations  is  also  found  in 
Guibas,  Ramshaw  and  Stolfi  (1983),  in  Fortime  (1985),  and  in  Edelsbrunner 
and  Welzl  (1985). 

The  problems  become  considerably  more  difficult  when  the  motion  of  Q 
may  also  involve  rotations.  Chazelle  (1983)  presents  an 
0(m^n^{m  +  n)  log  (m+n))  naive  algorithm  for  the  general  polygon 
containment  problem.  If  Q  is  a  line  segment  and  P  is  an  arbitrary  polygonal 
region,  then  the  general  motion  planning  problem  for  Q  can  be  solved  in 
^(n^log  n)  time  (cf.  Leven  and  Sharir,  1985a,  Sifrony  and  Sharir,  1985).  If 
C  is  a  convex  polygon  (and  P  again  is  an  arbitrary  polygonal  region)  then 
motion  of  Q  can  be  planned  in  time  0{mn\^{mn)  log  mn)  (Leven  and  Sharir, 
1985c,  Kedem  and  Sharir,  1986),  where  X,(/fe)  is  an  almost  linear  function  of 
k  yielding  the  maximal  number  of  connected  graph  portions  which  compose 
the  lower  envelope  of  k  continuous  fimctions,  each  pair  of  which  intersect  in 
at  most  s  points.  Chazelle  (1983)  has  shown  that  if  P  is  also  convex,  then  the 
existence  of  a  (translated  and  rotated)  copy  of  Q  inside  P  can  be  determined 


in  0{mn^)  time.  Another  recent  related  work  by  Yap  (1986)  involves 
planning  the  passage  of  an  arbitrary  simple  polygon  P  through  a  'door'  (an 
interval  opening  in  some  infinite  line  obstacle  /).  Such  a  motion  (which  can 
also  be  viewed  as  separation  of  P  from  another  polygonal  arc  Q,  chosen  so 
that  it  overlaps  /  in  a  sufficiently  large  interval,  and  its  endpoints  are  the 
endpoints  of  the  door)  can  be  planned  in  time  0{n^). 

As  stated  above,  in  this  paper  we  investigate  the  problem  of  separating 
Q  from  P  by  a  purely  translational  motion.  We  develop  an  algorithm  which 
runs  in  time  0{mn  a(mn)  log  m  log  n),  where  a(^)  is  the  inverse 
Ackermann's  function.  We  also  exhibit  an  example  in  which  Q  may  require 
n(mn)  translations  to  be  separated  from  P,  showing  that  in  the  worst  case 
our  algorithm  is  close  to  being  optimal. 

We  also  show  that  the  space  F  of  all  free  placements  of  Q  (whose 
combinatorial  complexity  is  always  at  most  0(m^n^))  can  have  Cl^m^n'^) 
connected  components  in  the  worst  case.  Thus  our  algorithm  is  superior  to 
any  motion-planning  algorithm  that  has  to  calculate  the  entire  space  F. 

Our  algorithm  has  also  the  following  properties: 

(1)  Given  a  final  desired  separated  position  of  Q ,  the  algorithm  can  produce 
the  shortest  separating  translational  motion  of  Q  from  its  given  position  to 
that  destination. 

(2)  Given  an  integer  k,  the  algorithm  can  determine  whether  Q  can  be 
separated  from  P  using  at  most  k  translations,  and,  if  so,  produce  such  a  "k- 
separating  motion". 
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(3)   The  algorithm  can  be  generalized  to  an  algorithm  for  planning  collision- 
free  translational  motion  of  Q  between  any  two  free  placements. 

2.   The  algorithm 

Chir  algorithm  is  based  on  the  following  well  known  observation  (see 
Lozano  Perez  and  Wesley,  1979).  Fix  a  reference  point  Z  inQ,  and  suppose 
without  Ibss  of  generality  that  at  the  given  placement  of  Q,  Z  lies  at  the 
origin.   Define 

K  =-  P  -  Q  =  {x-y  :x  €  P,y  €  Q} 
where  x—y  denotes  vector  difference  {K  is  known  as  the  Minkowski  or 

vector  difference  of  P  and  Q).    Qearly  a  placement  of  Q  (with  the  same 

given  orientation)  intersects  P  if  and  only  if  the  reference  point  Z  at  this 

placement  lies  in  K.  Thus  the  space  F  of  free  placements  of  Q  (that  is, 

placements  in  which  P  and  Q  are  disjoint)  is  conveniently  represented  as 

F  =  ii:*^  in  the  sense  that  each  x  €  F  corresponds  to  the  free  placement  of  Q 

in  which  Z  coincides  with  x. 

Thus  our  goal  is  reduced  to  that  of  analyzing  K".  More  specifically,  we 
wish  to  determine  whether  the  origin  (i.e.  the  given  placement  of  Z)  and  the 
point  at  infinity  lie  in  the  same  connected  component  of  K",  and,  if  so, 
calculate  a  (polygonal)  path  it  between  these  two  points  which  lies  entirely 
within  K''.  Note  that  each  straight  segment  in  tr  corresponds  to  a  single 
translation  of  |2.  so  that  the  number  of  segments  in  ir  is  the  nimiber  of 
translations  in  which  Q  can  be  separated  from  P;  also  the  length  of  ir 
between  O  and  some  point  sufficiently  far  from  P  is  equal  to  the  total 
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translational  distance  in  which  Q  is  moved  during  this  motion. 

The  approach  that  we  take  thus  aims  to  calculate  the  unbounded 
connected  component  €„  of  /sT".  We  first  establish  some  properties  of  K  and 
ofC„. 

Lemma  2.1:  (a)   AT  is  a  polygonal  region  having  at  most  0(m^n^)  comers, 
(b)  C«  is  a  simple  polygon. 

Proof:  (a)  We  repeat  here  well  known  arguments  (cf.  for  example  Guibas, 
Ramshaw  and  Stolfi,  1983).  Qearly  the  boundary  of  K  must  consist  of  vector 
differences  of  pairs  of  points  lying  respectively  on  the  boundaries  of  P  and  of 
Q.  Hence  P  —  Q  is  3.  polygonal  region,  each  of  whose  edges  has  the  form 
p  —  q,  where  either  p  is  an  edge  of  P  and  ^  is  a  vertex  of  2  or  P  is  a  vertex 
of  P  and  5  is  an  edge  of  Q.  Since  there  are  at  most  2mn  such  differences,  it 
follows  that  the  boundary  of  P  —  Q  is  contained  in  the  union  of  these  2mn 
segments.  Moreover,  each  comer  oi  P  —  Q  must  be  either  the  difference  of 
a  vertex  of  P  and  a  vertex  of  Q,  or  a  point  of  intersection  of  two  of  the 
above  segments.  Since  there  are  plainly  at  most  mn  comers  of  the  first  kind 
and  0{w}rp-)  comers  of  the  second  kind,  the  claim  follows. 

(b)  It  is  plain  that  P  —  Q  \s  connected,  so  that  the  boundary  of  any 
connected  component  of  {P  —  QY  must  be  connected.  □ 

Before  continuing  we  present  two  examples  which  help  to  calibrate  the 
worst  case  combinatorial  complexity  otP  -  Q  and  of  C,. 

Example  1:  This  example  shows  that  in  the  worst  case  P  —  Q  can  have 
^{w}n^)  connected  components  (and  thus  also  n(m^n^)  comers).    In  this 
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example,  as  illustrated  in  Fig.  2.1,  Q  consists  of  two  'combs'  each  having  m 
long  thin  'teeth'  so  that  their  'backbones'  are  perpendicular  to  one  another 
and  their  teeth  point  away  from  the  right  angle  formed  between  the 
backbones.  P  also  consists  of  a  pair  of  matching  aligned  combs,  each  having 
n  teeth,  which  now  point  into  the  right  angle  formed  between  the  backbones 
of  P.  The  separation  between  adjacent  teeth  of  F-is  taken  to  be  much  larger 
than  the  separation  between  adjacent  teeth  of  Q.  (Fig.  2.1  shows  a  skeletal 
representation  of  P  and  Q ;  by  slightly  thickening  these  skeletons  we  can  turn 
them  both  into  simple  polygons).  It  is  clear  that  if  the  length  of  the  teeth  and 
backbones  of  P  and  Q  and  the  separations  between  teeth  of  P  and  between 
teeth  of  Q  are  appropriately  chosen,  the  space  P—Q  of  free  positions  of  Q 
will  contain  n(m^/i^)  connected  components,  each  of  which  (except  the 
unbounded  component)  is  determined  by  choosing  a  pair  of  adjacent  vertical 
teeth  of  Q  and  placing  them  on  two  sides  of  a  vertical  tooth  of  P,  and  by 
choosing  similarly  a  placement  of  the  horizontal  teeth  of  Q  amidst  those  of  P. 

Example  2:  This  example  shows  that  in  the  worst  case  the  boundary  of  C„ 
can  have  Cl(mn)  comers.  In  this  example,  as  illustrated  in  Fig.  2.2(a),  P  has 
a  skeletal  representation  consisting  of  a  sequence  of  m'  rectangular  'rooms' 
lying  in  a  row  next  to  each  other,  such  that  each  pair  of  adjacent  rooms  are 
connected  by  a  small  'door'  in  their  common  wall,  and  such  that  the  last 
(rightmost)  room  has  also  a  similar  door  in  its  right  exterior  wall.  The 
second  polygon  Q  has  a  skeletal  representation  of  the  shape  of  a  zigzag  line 
consisting  of  n'  segments.  The  initial  placement  of  j2  is  in  the  'innermost' 
(leftmost)  room  of  P.  The  dimensions  of  P  and  Q  can  be  chosen  so  that  the 
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only  way  to  translate  Q  out  of  P  is  to  move  it  to  the  right  through  one  door 
at  a  time,  and  so  that  translation  of  Q  through  each  door  must  involve  n' 
distinct  translations  in  alternating  upward  and  downward  directions,  each 
pushing  a  different  segment  of  (the  skeleton  of)  Q  through  the  door.  Since  P 
has  m=0{ni')  sides  and  Q  has  n  =  0{n')  sides,  it  follows  that  in  this  example 
n(mn)  translations  may  be  required  to  separ^tisQ  from  P.  This  example  also 
shows  that  the  unbounded  component  C„  of  (.P,~QX  •''^^  consist  of  CL{mn) 
sides,  as  illustrated  in  Fig.  2.2(b).  We  are  indebted  to  Ryan  Hayward  for 
suggesting  this  example. 

We  next  show  that  even  though  the  entire  P-Q  may  have  Cl{m^n-^) 
corners,  the  boundary  bd(C„)  of  C„  contains  at  most  only  0{mn  a(mn)) 
comers,  a  bound  which  is  very  close  to  the  worst-case  lower  bound  provided 
by  Example  2. 

Theorem  2.2:  bd(C„)  has  at  most  0(mn  oc(mn))  comers. 

Proof:  As  observed  in  the  proof  of  Lemma  2.1(a),  bd{Cj)  is  contained  in  the 
union  of  2mn  segments,  each  of  which  is  either  a  difference  of  a  side  of  P 
and  a  vertex  of  j2  or  of  a  vertex  of  P  and  a  side  of  Q.  Enujnerate  these 
segments  as  ei,C2,  .  .  .  ,c,,  where  q=2mn.  It  is  clear  that  we  can  orient  each 
segment  e^  so  that  a  sufficiently  small  neighborhood  of  e,  lying  on  its  right 
side  is  disjoint  from  Co..  By  Lemma  2.1(b),  y  =  bd{C„)  is  a  simple  dosed 
polygonal  curve,  and  the  above  observation  implies  that  if  we  traverse  7  in  a 
clockwise  direction,  then  for  each  segment  e,,  every  portion  of  it  that  appears 
along  -y  is  traversed  in  the  direction  assigned  to  c<;  moreover,  the  (clockwise) 


order  in  which  these  portions  are  encountered  along  7  coincides  with  their 
order  along  e/  (we  omit  details  of  the  proof  of  these  rather  simple  topological 
facts). 

Consider  the  circular  sequence  of  straight  segments  71,72,  ...  ,7,  of  which  7 
is  composed  (starting  at  an  arbitrary  corner  of  7  and  arranged  in  clockwise 
order).  For  each  i^t  the-  segment  7,  is  a  portion  of  some  e„,;  let  U  denote 
the  (circular)  sequence  ui,U2.  — "•  \'u,:  For  each  segment  e^  consider  all 
appearances  of  its  index  a  in  U.  The  preceding  arguments  imply  that  there 
exist  two  appearances  w/(a),  «/(a)  of  a  in  U,  which  we  denote  respectively  as 
the  designated  first  and  the  designated  last  appearances  of  a  in  U,  such  that  all 
other  appearances  of  a  in  U  are  within  its  portion 
tA")  =  (  M/(a) ,«/(«)+!,  .  .  .  ,«/(a)- !,«/(«)  )•  (Regarding  these  notations,  recall 
that  U  is  a.  circular  sequence,  so  we  might  have/(a)>/(a),  in  which  case  U^'^^ 
consists  of  the  portion  Uff^^y  ...,«,  followed  by  the  portion  ui,  .  .  .  ,M;(a); 
note  also  that  one  might  have  f{a)=l{a),  in  which  case  a  appears  in  U  just 
once,  and  CA")  consists  of  the  single  element  "/(a).)  We  will  regard  each  U^"^ 
as  a  linear  sequence  whose  elements  are  ordered  by  their  circular  order  along 
U. 

The  sequence  U  has  the  following  properties: 

(1)  Ui  ^  {1,  .  .  .  ,q}  for  eadi  i^t. 

(2)  U/#«^+i  for  each  i^t  (where  /+ 1  is  computed  modulo  t). 

(3)  For  each  pair  a  i'  b  ^  {1,  .  .  .  ,q},  there  do  not  exist  five  indices 
p<i<j<k<q  in  CA")  (in  the  linear  order  induced  on  L^")  by  V)  such  that 
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Up=Uj  =  Uq  =  a  and  Ui  =  ui^  =  b. 

The  first  two  properties  are  obvious.  To  establish  the  third  property, 
suppose  to  the  contrary  that  such  a,  b,  and  p  ,ij,k,q  exist.  Let  us  denote  by 
7/ J  the  portion  of  y  traversed  in  clockwise  order  from  7,  to  yj,  excluding  y, 
and  yj.  Without  loss  of  generality,  suppose  that  no  element  of  U  between  Ut 
and  Ujfc  is  equal  to  b,  and  that  no  element  of  U  betv/een  Up  and  u,  or  between 
Ut  and  Ug  is  equal  to  a.  We  distinguish )^^k:^§SJ  .laa^.istf;, 

(i)  The  portion  e*  of  e^  between.^: 7^-  'and  7^  (excluding  these  two 
subintervals)  does  not  intersect  the  portion  e^  of  e^,  between  its  two 
subintervals  7/  and  7^  (again,  excluding  7,  and  7^).   Let  8  be  the  closed  curve 

8  =  7;,./   7/   e^   7jk   7;fc,^   e* 
(see  Fig.  2.3).    By  our  assumptions  and  the  fact  that  7  is  simple,  it  follows 

that  8  is  also  simple,  and  is  thus  a  (polygonal)  Jordan  curve. 

We  claim  that  the  interior  of  8  is  disjoint  from  Co.  Indeed,  suppose  the 
contrary  and  let  w  be  a  point  in  the  interior  of  8  which  also  belongs  to  C„. 
Then  one  can  find  a  path  ir  connecting  w  to  the  point  at  infinity,  lying 
entirely  within  the  interior  of  €„,  and  intersecting  8  only  finitely  many  times, 
so  that  each  of  these  intersections  is  transversal  and  takes  place  at  a  point  in 
the  relative  interior  of  some  edge  of  8.  But  each  such  intersection  x  is  either 
at  an  edge  of  7,  i.e.  a  point  on  the  boundary  of  C„,  or  is  a  point  on  e^  or  on 
ei, ,  in  which  case  it  must  contain  points  (lying  sufficiently  near  x  on  one  side 
of  it)  which  lie  in  C« ,  contradicting  in  both  cases  the  assumption  that  tt  lies 
entirely  within  the  interior  of  C,.. 
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Now  consider  7^1  which  is  a  portion  of  e* .  Since  points  lying  on  the  left 
side  of  7y  sufficiently  near  it  belong  to  C„,  it  follows  that  the  interior  of  5 
near  7/  lies  on  the  right  side  of  that  edge.  But  then  it  is  easily  checked  that 
both  sides  of  7^  (also  of  7^)  sufficiently  near  these  edges  lie  in  the  interior  of 
8 ,  which  is  impossible,  because  one  of  these  sides  contains  points  in  C„ .  This 
contradiction  completes  the  argument  in  case  (i) .  - 

(ii)  e*  and  e^  intersect.  'L&9W-d%iiiiit'\k^ii  point  of  intersection,  and  suppose 
without  loss  of  generality:that  x  4ies  between  7^  and  7^  (see  Fig.  2.4).  Let  8 
be  the  curve 

8  =  lp,i  7/  et,   e^  , 
where  e^  (resp.  F/,)  is  the  portion  of  e^  (resp.  of  e^,)  between  7^  and  x, 

excluding  7^  (resp.   between  7,  and  x,  excluding  7^).  Arguing  as  above,  it  is 

easy  to  show  that  8  is  a  Jordan  curve  whose  interior  is  disjoint  from  C.. 

However,  it  follows  from  the  structure  of  8  that  either  the  left  side  of  7/  or 

the  left  side  of  7^  (sufficiently  near  these  edges)  must  lie  in  the  interior  of  8 , 

which  is  impossible  because  these  sides  both  lie  in  Coo.  This  contradiction 

completes  the  proof  of  property  (3). 

We  next  transform  U  into  another  sequence  U*  as  follows.  Replace  each 
index  a  ^  q  for  which  U^"^  "wraps  around"  U  (i.e.  /(a)  >  l{a)),  by  two 
distinct  symbols  a' ,  a",  such  that  all  appearances  of  a  in  the  subsequence 
«/(a),  .  •  .  ,Ut  are  replaced  by  a' ,  and  all  appearances  of  a  in  ui,  .  .  .  ,«/(«)  are 
replaced  by  a" .  The  resulting  sequence  U*  has  the  same  length  t  as  U,  and  is 
composed  of  at  most  2q  symbols.  We  claim  that  U*  satisfies  the  modified 
properties 
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(2*)   u*¥=u*+i  for  each  i<t 

(3*)  For  each  pair  of  symbols  a  ¥=  b,  there  do  not  exist  five  indices 
/i<j2<i3<J4<i5  in  U*  such  that  u^*  =  u^*  =  U/*  =  a  and  ui*  =  ut*=b. 

Indeed,  concerning  (3*),  this  property  is  a  direct  consequence  of  (3)  if  a  is 
one  of  the  original  unsplit  indices  of  U.  If,  say,  a  =  c'  for  some  original 
index  c,  then  all  indices  j'l,  .  ,  .  ,;5  must  also  belong  to  the  portion  of  lA'^^ 
between  the  indices  /(c)  and  t,  so  agnitrthe  claimi "fellows  from  property  (3) 
of  U.  Similar  arguments  apply  if  a  ~-  c"  for  some,  ind^x  c. 

U*  is  thus  a  (2q ,3)-DavenportSchinzel  sequence  in  the  terminology  of  Hart 
and  Sharir  (1984),  and  by  the  results  of  that  paper,  the  length  of  U*  (and  of 
U)  is  at  most  0(2q  a{2q))  =  0(mn  a(mn)).  n 

Remark:  Using  similar  arguments  to  those  in  the  proof  of  Theorem  2.2,  one 
can  obtain  the  following  generalization:  Let  e^,  .  .  .  ,e„  be  any  (possibly 
intersecting)  n  straight  segments  in  the  plane.  Then  the  boundary  of  the 
unbounded  (in  fact,  of  any)  component  of  the  complement  of  the  union  of 
these  segments  consists  of  at  most  0(n  a(n))  segments  (which  are  portions  of 
the  segments  e,). 

Efficient  calculation  of  bd{C^) 

We  next  present  an  efficient  algorithm  for  the  calculation  of  y  =  bd{C,o)- 
First  obtain  a  hierarchical  decomposition  of  (2  as  in  Chazelle  (1982). 
Specifically,  we  first  obtain  a  triangulation  T  of  Q.  Then  we  cut  Q  along  one 
of  the  diagonals  of  T  so  as  to  divide  it  into  two  subpolygons  Qi,  Qi  each  of 
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which  contains  at  least  some  fixed  fraction  of  the  sides  of  Q ,  and  continue  to 
cut  j2i.  Qi  recursively  in  the  same  manner.  Chazelle  (1982)  has  shown  that 
such  a  decomposition  is  always  possible;  Guibas  et  al.  (1986)  present  a  linear 
time  algorithm  for  the  calculation  of  such  a  decomposition. 

Our  algorithm  then  applies  the  following  divide  and  conquer  approach: 
Let  Q=QiUQ2  be  the  d©S»njposition  of  Q  as  described  above.  Calculate 
recursively  the  boundaries  ;iy!i^  sfi^o^'Ahe  unbounded  connected  component  of 
(P—QiY,  (P  —  QiY  respectively.'  Then  merge  yi,  72  to  obtain  the  desired 
boundary  7  =  bd(C„)  of  the  unbounded  component  of  {P—QY.  The 
merging  of  71  and  72  is  perfonned  as  follows.  Since  71  and  72  are  both 
simple  polygons,  we  can  use  the  technique  of  Chazelle  and  Guibas  (1985) 
(see  also  Guibas  et  el.  1986)  to  preprocess  each  of  these  polygons  into  a 
data-structure  which  supports  efficient  response  to  ray  shooting  queries, 
where  each  such  query  asks  for  determination  of  the  first  point  on  71  (resp. 
on  72)  hit  by  a  ray  emerging  from  some  specified  point  X  in  a  specified 
direction  u.  As  shown  in  Chazelle  and  Guibas  (1985)  and  in  Guibas  et  al. 
(1986)  this  shooting  problem  can  be  solved  in  0(t)  preprocessing  time  (and 
storage)  and  0(log  t)  query  time,  where  t  is  the  number  of  sides  of  71  (resp. 
of  yi)- 

Having  preprocessed  71  and  72  in  this  manner,  we  next  find  a  starting 
point  Xi  lying  on  one  of  these  curves  and  being  exterior  to  the  other  (e.g. 
one  can  take  X^  to  be  the  leftmost  vertex  among  all  vertices  of  71,  72).  We 
then  begin  to  trace  the  desired  boundary  7,  which  is  also  easily  seen  to  be  the 
boundary  of  the  imbounded  connected  component  of  (7iU72)*^,  from  Xi  in  a 
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clockwise  direction.  Suppose  we  have  progressed  along  y  from  Xi  up  to 
some  comer  A,  of  7.  If  AT,  is  a  comer  of  7^  or  of  72.  say  for  definiteness  of 
7i,  we  take  the  next  edge  of  y  to  be  traced  to  be  the  edge  of  7^  incident  to  Xi 
and  pointing  from  it  in  clockwise  direction  along  71.  Suppose  on  the  other 
hand  that  Xi  is  a  point  of  intersection  of  an  edge  e^  of  71  with  an  edge  62°^ 
72-  Suppose  further  without  loss  of  generality  that  -oi^r  tracing  of  7  has 
reached  Xi  along  62',  then  the  next  edge  e  ^-i-v^o^l^e  Jtaced  is  taken  to  be  ei 
(which  has  to  be  traced  from  Xi  in  clockwise  direction  along  71). 

In  either  case  we  now  have  an  edge  e  of  71  (and  of  7)  which  we  want  to 
trace  from  the  point  X,  lying  on  it  until  the  next  comer  X/+i  of  7.  To  find 
Xi+i  we  perform  a  ray  shooting  query  to  find  the  first  point  Z  on  72  hit  by 
the  ray  emerging  from  AT,  in  the  direction  of  (the  appropriate  portion  of)  e. 
If  Z  lies  on  e,  we  put  ^,+1  :=  Z;  otherwise  we  take  X,+i  to  be  the  appropriate 
endpoint  of  e.  This  tracing  process  is  repeated  until  we  trace  the  complete 
boundary  7  back  to  Xi. 

The  complexity  of  this  merging  procedure  can  be  estimated  as  follows. 
Let  n,  be  the  number  of  sides  of  Qi  for  i  =  l,2.  By  Theorem  2.2,  7,  has  at 
most  0(mnia{mni))  vertices,  so  that  the  preprocessing  of  71  and  72  for  the 
shooting  queries  is  accomplished  in  overall  time 

0(mnia(mni))  +  0(mn2a(mn2))  =  0{mn  a(mn))  . 
The  starting  point  Xi  can  also  be  calculated  in  0{mn  a{mn))  time.    The 

tracing  of  7  consists  of  repeated  applications  of  ray  shooting  queries,  one  for 

each  comer  of  7.   Since  7  has  at  most  0{mn  a(mn))  comers,  it  follows  that 

the  complexity  of  tracing  7  is  0{mn  a(mn)  log  mn). 
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Let  T{m,n)  denote  the  maximal  time  required  to  calculate  the  boundary 
of  the  unbounded  component  of  the  complement  of  P—Q,  where  P  and  Q  are 
simple  polygons  having  m  and  n  sides  respectively.  Then  we  have  the 
following  recurrence: 

T{m,n)  ^  T{m,n{)  +  T{m,n2)  +  0{mn  a{mn)  log  mn) 
where  n>3  and  where  both" wi-aiid  ^2  are  ^^  (cf.  Chazelle,  1982,  Guibas  et 

al.,  1986).  Thisf^c^rfetiik^fti^'li^  ^  ^§^^  '"''-' 

T(m,n)  =  0(nT{m,3))  +  0{mn  a(mn)  log  mn  log  n)  . 
To  calculate  r(m,3)  we  make  use  of  the  fact  that  in  this  case  Q  is  just  a 

triangle,    and    is    therefore    convex.     Thus    calculation    of   P—Q    can    be 

accomplished  by  the  generalized  Voronoi  diagram  approach  of  Leven  and 

Sharir  (1985b)  in  time  0{m  log  m).   It  follows  that 

T(m,n)=0(mn  [  log  m  +  a(mn)  log  mn  log  n  ])  , 
and,  assvuning  without  loss  of  generality  m  ^  n,  we  obtain 

T{m,n)  =  0{mn  a(mn)  log  m  log  n)  . 

Remarks  and  Open  Problems 

(1)  The  bound  0{mn  a(mn))  on  the  size  of  bd{C»)  as  given  by  Theorem  2.2 
is  not  known  to  be  tight  in  the  worst  case;  the  best  matching  lower  boimd  is 
that  given  in  Example  2.  Moreover,  for  n  =  3  or,  more  generally,  for  a 
convex  Q,  it  is  known  that  the  entire  P—Q  contains  0{mn)  corners  of  which 
only  0{m)  comers  are  formed  by  intersection  of  edges  (as  in  the  proof  of 
Lemma  2.1(a);  cf.  Kedem  and  Sharir,  1985,  Kedem  et  al.,  1986,  Leven  and 
Sharir,  1985b).   Although  for  all  practical  purposes  the  bound  0{mn  a(mn)) 
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can  be  considered  to  be  the  same  as  0{mn),  it  is  still  interesting  from  a 
theoretical  point  of  view  to  improve  it,  at  least  in  some  special  cases,  e.g. 
when  Q  can  be  decomposed  as  the  union  of  a  small  number  of  openly  disjoint 
convex  subpolygons.  Alternatively,  can  one  show  that  in  the  worst  case 
bd{C„)  can  indeed  consist  of  Q,{mn  a(mn))  corners? 

(2)  The  merging  procedure  described  auove  deserves  sbrne  comments.  First 
of  all,  it  is  significant  that  this  probediif^(^'d([:^' mjf  c^ciulate  all  intersections 
between  y-i  and  72  (of  which  there  ruight  coiiceivably  be  n(m^n^)).  Note  also 
that  if  7i,  -y2  or  7  were  not  connected,  then  the  complexity  of  our  procedure 
would  deteriorate,  either  because  we  would  have  to  spend  more  time  in 
locating  a  starting  point  on  each  component  of  7,  or  because  the  shooting 
queries  would  require  more  time,  because  the  technique  of  Chazelle  and 
Guibas  (1985)  applies  only  to  simple  polygons.  Finally,  it  is  still  an  open 
problem  whether  7  can  be  calculated  in  time  linear  in  the  nimiber  of  sides  of 
7i,  72  and  7. 

(3)  Theorem  2.2,  and  the  calculation  of  bd(C„)  can  be  generalized  to  apply 
to  any  connected  component  of  {P-QY.  In  particular,  we  can  apply  them  to 
the  connected  component  C  of  {P—QY  which  contains  the  given  position  of 
(the  reference  point  Z  on)  Q.  Thus  in  time  0{mn  a(mn)  log  m  log  n)  we  can 
calculate  the  space  of  all  free  placements  of  Q  reachable  from  its  given 
position  by  collision-free  translational  motion. 

3.  Calculating  Motions  Separating  Q  from  P 
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In  this  section  we  comment  on  the  actual  calculation  of  the  desired 
translational  motion  separating  Q  from  P,  using  the  data  structures  calculated 
in  the  previous  section.  Let  y  =  bd{Coo)  be  the  (simple  polygonal)  boundary 
of  the  unbounded  connected  component  C„  of  (P  —  QY,  and  let 
t=0{mn  a(mn))  be  the  number  of  sides  of  7.  To  determine  whether  Q  can 
be  separated  from  P  b^  a  p^ijcelyytj^slational  motion  from  a  given  placement, 
we  simply  have  to^  .t^  y^^tj^jtl^e j^^jsjti^  O  of  the  reference  point  on  Q  at 
that  placement  Hes  ija  ^p,..^This^ea|]i^',b,6„ti3vially  done  in  0(t)  time.  If  O  does 
indeed  lie  in  C„,  we  can  calculate  an  'optimal'  translational  separating  motion 
of  Q  in  the  sense  of  one  of  the  following  two  approaches: 

(i)  Calculate  a  Euclidean  shortest  separating  motion:  Let  7*  be  the  convex 
hull  of  7.  We  will  follow  the  reasonable  convention  that  Q  and  P  are 
considered  to  be  fully  separated  whenever  O  lies  on  or  outside  7*.  Again  we 
can  test  whether  O  lies  on  or  outside  7*  in  0(t)  time.  If  so,  no  separating 
motion  is  required.  Otherwise  O  must  lie  in  one  of  the  'pockets'  enclosed 
between  7  and  7*,  where  each  such  pocket  i?  is  a  simple  polygonal  region 
bounded  by  an  edge  e^  of  7*  (the  'lid'  of  R)  not  belonging  to  7  and  by  a 
contiguous  portion  of  7  (see  Fig.  3.1).  In  this  case  our  goal  is  to  calculate  the 
shortest  path  within  R  from  C?  to  a  point  in  e^.  This  can  be  done  in  0(tii) 
time,  where  t^  is  the  number  of  sides  of  R,  using  the  shortest  path  algorithm 
of  Guibas  et  al.  (1986).  This  algorithm  calculates  the  shortest  paths  from  O 
to  all  vertices  of  R,  and  produces  a  partitioning  of  R  into  triangular  regions 
such  that  for  eadi  of  these  triangles  A  there  exists  a  vertex  v  of  R  such  that 
the  terminal  segments  of  the  shortest  paths  from  O  to  all  the  points  in  A  all 
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emerge  from  v.  Hence  by  examining  each  triangle  in  this  partitioning  which 
intersects  ej^,  it  is  straightforward  to  calculate  the  desired  shortest  path  to  e^f 
in  0(//j)  time. 

(11)  Calculate  a  separating  motion  consisting  of  the  smallest  number  of 
translations:  This  is  also  considered  to  be  a  useful  criterion  for  efficiency  of 
the  separating  motion  (cf.  Toussaint,  1.9^5}i  gWfetFHlsay  that  P  and  Q  are  k- 
sep arable  by  translations  if  ♦.h^yoiC^Wq  beti^.^^^1^  t)y  a  sequence  of  k 
translations,  but  not  by  any  sequence,  of- fewOT-translations;  such  a  sequence 
of  translations  will  be  called  a  k-separating  motion  (of  Q  from  P).  To  find 
such  a  ^-separating  motion,  we  first  test,  as  in  (i)  above,  whether  O  lies 
outside  7*,  in  which  case  k=Q  and  no  motion  is  required.  If  O  lies  inside  7*, 
let  R  be  the  pocket  containing  O  (as  in  (i))  and  let  e^  be  its  lid.  Our  task  is 
then  to  find  a  polygonal  path  within  R  from  O  to  e^  consisting  of  the  fewest 
possible  number  of  edges.  This  problem  has  been  recently  studied  by  Suri 
(1986).  To  describe  his  results,  let  us  partition  R  into  a  collection  Vi{en), 
1  =  1,2,  ,  .  .  ,  of  polygonal  regions  defined  as  follows.  Vi{en)  consists  of  all 
points  in  R  directly  visible  (within  R)  from  some  point  on  e^.  Inductively, 
^/+i(^j?)  consists  of  all  points  in  R  which  are  visible  from  some  point  in 
Vi{eji)  and  are  not  contained  in  U  VAcr).   Suri  shows  that  this  partitioning  of 

R  can  be  calculated  in  time  0{tn).  Given  this  partitioning,  we  find  (also  in 
0{tj()  time)  the  region  V)t(«/j)  containing  O.  Then  clearly  P  and  Q  are  k- 
separable  from  one  another,  and  the  data  structures  produced  by  Suri's 
algorithm  enable  us  to  calculate  a  iS: -separating  motion  of  Q  also  in  0(r/e) 
time. 
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In  summary,  we  have  shown 

Theorem  3.1:  After  calculating  7,  one  can  determine  in  additional  time  O(r^) 
(where  t^  is  the  number  of  comers  of  7)  whether  Q  can  be  separated  from  P 
by  translations,  and,  if  so,  also  calculate  such  a  separating  motio  either 
having  a  minimum  length  or  a  minimum  number  of  links. 

Remark:  The  twoapproadi'es  just  described  can 'be  easily  modified  so  that 
they  first  perform: '4  ^i^fSQ^s&i^  pfeSse,. Which  depends  only  on  the  shape  of 
P  and  Q,  but  not  dn  theif-jilFeseht  placements,  and  then,  given  specific 
placements  of  P  and  Q ,  determine  quickly  whether  translational  separability 
of  Q  from  P  from  these  placements  is  possible,  and,  if  so,  also  calculate  the 
shortest  Euclidean  length  of  such  a  separating  motion,  or,  alternatively,  the 
smallest  number  of  links  in  such  a  motion.  Using  the  techniques  of  Guibas  et 
al.  (1986)  and  of  Suri  (1986),  such  a  preprocessing  can  be  done  in  0{t)  time 
and  space,  and  each  actual  separability  query  can  be  answered  in  0(log  t) 
time. 
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Fig.  1.1.  An  instance  of  tlie  polygon  separability  problem. 


Fig.  2.1.  Example  1. 


-23- 


Fig.  2.2(a).  Example  2. 


Fig.  2.2(b).  (P  -  QY  =  C«  in  Example  2. 
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Fig.  2.3.  Case  (i)  in  the  proof  of  Theorem  2.2,  property  (3). 


Fig.  2.4.  Case  (ii)  in  the  proof  of  Theorem  2.2,  property  (3). 
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Fig.  3.1.  The  curves  7,7*,  and  a  shortest  separating 
motion  of  Q  from  a  position  within  some  pocket  R. 


,  book  roriy  be  kept 

FOURTEEN 

.,  ,lav  the  book  is  kept  wertim... 
ne  «ill  be  charRed  for  c:>ch  >lay  the 


DAYS 


Til 


UORD   1*Z 


NYU  COMPSCI  TR-215 
Pollack,  Richard 
Separating  two  simple 
polygons  by  a  sequence  of 
translations       c.2 


NYU  COMPSCI  TR-215 
Pollack,  Richard 
Separating  two  simple 
polygons  by  a  sequence  of 
translations       c.2 


LIBRARY 

N.Y.U.  Courant  Institute  of 

Mathematical  Sciences 

251  Mercer  St. 
New  York,  N.  Y.    10012 


^ 


m 


■^ 


, 


i 


